A systematic ab-initio study of curvature effects in carbon nanotubes 
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We investigate curvature effects on geometric parameters, energetics and eiectronic structure of 
zigzag nanotubes with fuliy optimized geometries from first-principle calculations. The calculated 
curvature energies, which are inversely proportional to the square of radius, are in good agreement 
with the classical elasticity theory. The variation of the band gap with radius is found to differ from 
simple rules based on the zone folded graphene bands. Large discrepancies between tight binding 
and first principles calculations of the band gap values of small nanotubes are discussed in detail. 
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I. INTRODUCTION 

Single wall carbon nanotubes (SWNTs) are basically 
rolled graphite sheets, which are characterized by twa 
integers in, m) defining the rolling vector of graphitc.EJ 
Therefore, electronic properties of SWNTs, at first or- 
der, can be deduced from that of graphene by map- 
ping thcuhajad- structure of 2D hexagonal lattice on a 
cylinder BBoQiJ Such analysis indicates that the (n, n) 
armchair nanotubes are always metal and exhibit one di- 
mensional quantum conduction!! The (n, 0) zigzag nan- 
otubes are generally semiconductor and only are metal 
if n is an iateger multiple of three. However, recent 
experimental] indicate much more complicated structural 
dependence of the band gap and electronic properties of 
SWNTs. The semiconducting behavior of SWNTs has 
been of particular interest, since the electronic properties 
can be controlled by doping or implementi ng , diffc ct^in a 
nanotube-based optoelectronic devices. BBll^Olld't 3 ]^ It 
is therefore desirable to have a good understanding of 
electronic and structural properties of SWNTs and the 
interrelations between them. 

Band calculations of SWNTs were initially performed 
by using a one-band 7r-orbita l tightjbin ding modelfl Sub- 
sequently, experimental dataMliSliZEa on the band gaps 
were extrapolated to confirm the inverse proportional- 
ity with the radknR of the nanotube.u Later, first prin- 
ciples calculationlla within Local Density Approximation 
(LDA) showed that the a*-ir* hybridization becomes sig- 
nificant at small R (or at high curvature). Such an effect 
were not revealed by thepWpQstMtal tight-binding bands. 
Recent analytical studiesE30'E3 showed the importance 
of curvature effects in carbon nanotubes. Nonetheless, 
band calculations performed by using different methods 
have been at variance on the values of the band gap. 
While recent studies predict interesting effects, such_as 
strongly local curvature dependent chemical reactivityEJ, 
an extensive theoretical analysis of the curvature effects 
on geometric and electronic structure has not been car- 
ried out so far. 

In this paper, we present a systematic ab-initio analysis 




FIG. 1: A schematic side view of a zigzag SWNT, indicating 
two types of C-C bonds and C-C-C bond angles. These are 
labeled as di, di, Q\ and 62- Radius dependence of these vari- 
ables are important in tight-binding description of SWNTs as 
discussed in the text. 



of the band structure of zigzag SWNTs showing interest- 
ing curvature effects. Our analysis includes a large num- 
ber of zigzag SWNTs with n ranging from 4 to 15. The 
fully optimized structural and electronic properties of 
SWNTs are obtained from extensive first-principle calctu 
lations within the Generalized Gradient ApproximatiouE 3 ] 
(GGA) by using pseudopotential planewave mcthocO. 
We used plane waves up_to an energy of 500 eV and ul- 
trasoft pseudopotcntialaHj. The calculated total energies 
converged within 0.5 meV/atom. More details about the 
calculations can be found in Ref's[E6 27 . 



II. GEOMETRIC STRUCTURE 

First, we discuss effects of curvature on structural pa- 
rameters such as bond lengths and angles. Figure [l] shows 
a schematic side view of a zigzag SWNT which indicates 
two types of C-C bonds and C-C-C bond angles, re- 
spectively. The curvature dependence of the fully opti- 
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FIG. 2: (a) Normalized bond lengths (di/do and d,2/do) versus 
the tube radius R. (d = 1.41 A), (b) The bond angles (#i 
and 62) versus R. (c) The curvature energy, E c per carbon 
atom with respect to graphene as a function of tube radius. 
The solid lines are the fit to the data as 1 / R 2 . 



mized structural parameters of zigzag SWNTs are sum- 
marized in Fig. |2| The variation of the normalized bond 
lengths (i.e. dc-c/do where do is the optimized C-C 
bond length in graphene) and the bond angles with tube 
radius R (or n) are shown in Fig. ^a and b, respectively. 
Both the bond lengths and the bond angles display a 
monotonic variation and approach the graphene values 
as the radius increases. As pointed out earlier for the 
armchair SWNTsEa, the curvature effects, however, be- 
come significant at small radii. The zigzag bond angle 
(81) decreases with decreasing radius. It is about 12° less 
than 120°, namely the bond angle between sp 2 bonds of 
the graphene, for the (4, 0) SWNT, the smallest tube we 
studied. The length of the corresponding zigzag bonds 
(^2), on the other hand, increases with decreasing R. On 
the other hand, the length of the parallel bond (di ) de- 
creases to a lesser extent with decreasing R, and the angle 
involving this bond (#2) is almost constant. 

An internal strain is implemented upon the formation 
of tubular structure from the graphene sheet. The associ- 
ated strain energy, which is specified as the curvature en- 
ergy, E c , is calculated as the difference of total energy per 
carbon atom between the bare SWNT and the graphene 



(i.e. E c = E t .swnt ~ E Ttgraphene ) for 4 < n < 15. 
The calculated curvature energies are shown in Fig. ^c. 
As expected E c is positive and increases with increasing 
curvature. Consequently, the binding (or cohesive) en- 
ergy of carbon atom in a SWNT decrease with increas- 
ing curvature. We note that in the classical theory of 
elasticity thfi-cuxvature energy is given by the following 
sionBBEl 



expressior 



Yh 3 n 
~~24~R? 



a 

w 



(1) 



Here Y is the Young's modulus, h is the thickness of the 
tube, and f2 is the atomic volume. Interestingly, the ab- 
initio curvature energies yield a perfect fit to the relation 
a/R 2 as seen in Fig. |^c. This situation suggests that the 
classical theory of elasticity can be used to deduce the 
elastic properties of SWNTs. In this fit a is found to be 
2.14 eV A 2 /atom, wherefrom Y can be calculated with 
an appropriate choice of h. 



III. ELECTRONIC STRUCTURE 

An overall behavior of the electronic band structures 
of SWNTs has iteen revealed from zone folding of the 
graphene bandsfltra Accordingly, all (n, 0) zigzag SWNT 
were predicted to be metallic when n is multiple of 3, 
since the double degenerate tt and tt* states, which over- 
lap at the if -point of the hexagonal Brillouin ff»e (BZ) 
of graphene folds to the T point of the tubea'Q. This 
simple picture provides a qualitative understanding, but 
fails to describe some important features, in particular 
for small radius or metallic nanotubes. This is clearly 
shown in Table where the band gaps calculated in the 
present study are summarized and compared with results 
obtained from other methods in the literature. For ex- 
ample, our calculations result in small but non-zero en- 
ergy band gaps of 93, 78 and 28 meV for (9,0), (12,0) 
and (15,0) SWNTs, respectively (see Table 1). Recently, 
these gaps are measured by-. Scanning Tunneling Spec- 
troscopy (STS) experiments?! as 80, 42 and 29 meV, in 
the same order. The biggest discrepancy noted in Table | 
is between the tight-binding and the first-principles val- 
ues of the gaps for small radius tubes such as (7, 0). These 
results indicate that curvature effects are important and 
the simple zone folding picture has to be improved. More- 
over, the analysis of the LEI A bands of the (6, 0) SWNT 
calculated by Blase et a&B brought another important 
effect of the curvature. The antibonding singlet tt* and 
a* states mix and repel each other in curved graphene. 
As a result, the purely tt* state of planar graphene is 
lowered with increasing curvature. For zigzag SWNTs, 
the energy of this singlet 7r*-state is shifted downwards 
with decreasing R (or increasing ciixvature). Here, we ex- 
tended the analysis of Blase et ao to the (n, 0) SWNTs 
with 4 < n < 15 by performing GGA calculations. 

In Fig. |a, we show the double degenerate 7r-states 
(which are the valence band edge at the T-point), the 
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TABLE I: Band gap, E g , as a function of radius R of (n,0) zigzag nanotubes. M denotes the metallic state. Present results for 
E g were obtained within GGA. First row of Ref. [19] is LDA results while all the rest are tight-binding (TB) results. Two rows 
of Ref. B3] are for two different TB parametrization. 
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FIG. 3: (a) Energies of the double degenerate 7r-states (VB), 
the double degenerate 7r*-states (CB) and the singlet 7r*-state 
as a function of nanotube radius. Each data point corresponds 
to n ranging from 4 to 15 consecutively, (b) The calculated 
band gaps as a function of the tube radius shown by filled 
symbols. Solid (dashed) lines are the plots of Eq. §r(jE*h-|2|) 
The experimental data are shown by open diamondsu't 



energy of the singlet 7r*-state at the T -point of the BZ 
is above the doubly degenerate it* states {i.e. bottom of 
the conduction band), while it falls between the valence 
and conduction band edges for n = 7, 8, and eventu- 
ally dips even below the double degenerate valence band 
7r-states for the zigzag SWNT with radius less than 2.7 
A (i.e. n < 7). Therefore, all the zigzag tubes with ra- 
dius less than 2.7 A are metallic. For n = 7, 8, the edge 
of the conduction band is made by the singlet 7r*-state, 
but not by the double degenerate 7r*-state. The band 
gap derived from the zone folding scheme is reduced by 
the shift of this singlet tt* -state as a result of curvature 
induced a* — 7r* mixing. This explains why the tight 
binding calculations predict band gaps around 1 eV for 
n = 7, 8 tubes while the self-consistent calculations pre- 
dict much smaller value. 

Another issue we next address is the variation of the 
band gap, E g , as a function of tube radius. Based on 
7r-orbital tight binding model, it was proposedQ that E g 
behaves as 



tp dt) 
9 ^ l0 ~R~' 



(2) 



which is independent from helicity. Within the simple 
7r-orbital tight binding model, j is taken to be equal to 
the hopping matrix element V pp7! . (do is the bond length 
in graphene). However, as seen in Fig. |b, the band gap 
displays a rather oscillatory behavior up to radius 6.0 
A. The relation given in Eq. |2| was obtained by a sec- 
ond order Taylor expansion of .one-electron eigenvalues 
of 7r-orbital tight binding modelEl around the iC-point of 
the BZ, and hence it fails to represent the effect of the 
helicity. By extending the Taylor expansion to the next 
higher order, Yorikawa and MuramatsuE3'li3 included an- 
other term in the empirical expression of the band gap 
variation, 



double degenerate 7r*-states (which become the conduc- 
tion band edge at T for large R), and the singlet 7r*-state 
(which is in the conduction band for large R). As seen, 
the shift of the singlet 7r*-state is curvature dependent, 
and below a certain radius determines the band gap. For 
tubes with radius greater than 3.3 A (i.e. n > 8), the 



E a = V t 



do . 



PPT J£ [1 



-l)Vos(30) 



do, 
R l 



(3) 



which depends on the chiral angle, 9, as well as an in- 
dex p. Here 7 is a constant and the index p is defined 
as the integer from k = n — 2m = 3q + p. The fac- 
tor (— l) p comes from the fact that the allowed k is 
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nearest to either the K- or K - point of the hexago- 
nal Brillouin zone. For zigzag nanotubes studied here, 
the chiral angle is zero, so the second term just gives 
Br 2 dependence as ±^iV pp7T (dQ / R) 2 . Hence, the solid 
lines in Fig. ^|b are fits to the empirical expression. 
Eg = VppTrdo/R ± Vppn'ydl/R 2 , obtained from Eq. || 
for 6 =0 by using the parameters Vp Vlx = 2.53 eV-and. 
7 = 0.43. The experimental data obtained by STSE]'E!l 
are shown by open diamonds in the same figure. The 
agreement between our calculations and the experimental 
data is very good considering the fact that there might 
be some uncertainties in identifying the nanotube (i.e. 
assignment of (n, m) indices) in the experiment. The fit 
of this data to the empirical expression given by Eq. || 
are also presented by a dashed line for comparison. 

The situation displayed in Fig. ^ indicates that the 
variation of the band gap with the radius is not sim- 
ply 1/R, but additional terms incorporating the chirality 
dependence are required. Most importantly, the mixing 
of the singlet 7r*-state with the the singlet c*-state due 
to the curvature, and its shift towards the valence band 
with increasing curvature is not included in neither the 



7r-orbital tight binding model, nor the empirical relations 
expressed by Eq. || and 0. This behavior of the singlet 
7r*-states is of particular importance for the applied ra- 
dial deformation that modifi es the curvature and in turn 
induces metallizatior£3c3tJ. 

In conclusion, we investigated structural and electronic 
properties that result from the tubular nature of the 
SWNTs. The first-principles total energy calculations 
indicated that significant amount of strain energy is im- 
plemented in a SWNT when the radius is small. How- 
ever, the elastic properties can be still described by the 
classical theory of elasticity. We showed how the singlet 
7r*-state in the conduction band of a zigzag tube moves 
and eventually enters in the band gap between the dou- 
bly degenerate 7r* -conduction and 7r-valence bands. As 
a result, the energy band structure and the variation of 
the gap with radius (or n) differs from what one derived 
from the zone folded band structure of graphene based 
on the simple tight binding calculations. 

Acknowledgments This work was partially sup- 
ported by the NSF under Grant No. INT01-15021 and 
TUBITAK under Grant No. TBAG-U/13(101T010). 



1 M.S. Dresselhaus, G. Dresselhaus and P.C. Eklund, Sci- 
ence of Fullerenes and Carbon Nanotubes; Academic Press, 
San Diego (1996); R. Saito, G. Dresselhaus and M. S. Dres- 
selhaus, Physical Properties of Carbon Nanotubes; Impe- 
rial College Press, London (1998). 

2 N. Hamada, S. Sawada and A. Oshiyama, Phys. Rev. Lett. 
68, 1579 (1992). 

3 M.S. Dresselhaus, G. Dresselhaus and R. Saito, Phys. Rev. 
B 45, 6234 (1992). 

4 J.W. Mintmire, B.I. Dunlap and C.T. White, Phys. Rev. 
Lett. 68, 631 (1992). 

5 C.T. White, D.H. Robertson and J.W. Mintmire, Phys. 
Rev. B 47, 5485 (1993). 

6 S. Frank, P. Poncharal, Z. L. Wang and W. A. Heer, Sci- 
ence 280, 1744 (1998); a recent review: S. Ciraci, A. 
Buldum and I. Batra, J. Phys: Condens. Matter 13, 537 
(2001). 

7 M. Ouyang, J. Huang, C.L. Cheung and CM. Lieber, Sci- 
ence 292, 702 (2001). 

8 L. Chico, M.P. Lopez Sancho and M.C. Munoz Phys. Rev. 
Lett. 81, 1278 (1998). 

9 P.G. Collins, A. Zettl, H. Bando, A. Thess and R.E. Smal- 
ley Science 278, 100 (1997). 

10 M. Bockrath, D.H. Cobden, P.L. McEuen, N.G Chopra, 
A. Zettl, A. Thess and R.E. Smalley, Science 275, 1922 
(1997). 

11 A. Bezryadin, A.R.M. Verschueren, S.J. Tans, and C. 
Dekker, Phys. Rev. Lett. 80, 4036 (1998). 

12 g. Kdig, S. Ciraci, O. Giilseren and T. Yildirim, Phys. Rev. 
B 62, 16345 (2000). 

13 H.S. Sim, C.J.Park and K.J. Chang, Phys. Rev. B 63, 
073402 (2001). 

14 O. Giilseren, T. Yildirim and S. Ciraci, Phys. Rev. Lett. 
87, 116802 (2001). 

15 J.W.G. Wildoer, L.C. Venema, A.G Rinzler, R.E. Smalley 
and C. Dekker, Nature 391, 59 (1998). 



16 L.C. Venema, J.W. Janssen, M.R. Buitelaar, J.W.G. 
Wildoer, S.G Lemay, L.P. Kouwenhoven and C. Dekker, 
Phys. Rev. B 62, 5238 (2000). 

17 T.W. Odom, J. Huang, P. Kim and CM. Lieber, Nature 
391, 62 (1998). 

18 T.W. Odom, J. Huang, P. Kim and CM. Lieber, J. Phys. 
Chem. B 104, 2794 (2000). 

19 X. Blase, L.X. Benedict, E.L. Shirley and S.G. Louie, Phys. 
Rev. Lett. 72, 1878 (1994). 

20 C.L. Kane and E.J. Mele, Phys. Rev. Lett. 78, 1932 (1997). 

21 A. Kleiner and S. Eggert, Phys. Rev. B 63, 073408 (2001). 

22 A. Kleiner and S. Eggert, Phys. Rev. B 64, 113402 (2001). 

23 J. P. Perdew and Y. Wang, Phys. Rev. B 46, 6671 (1992). 

24 M.C. Payne, M.P. Teter, D.C. Allen, T.A. Arias and J. D. 
Joannopoulos, Rev. Mod. Phys. 64, 1045 (1992). 

25 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

26 T. Yildirim, O. Giilseren, C. Kihc and S. Ciraci, Phys. 
Rev. B 62, 12648 (2000); T. Yildirim, O. Giilseren and S. 
Ciraci, Phys. Rev. B 64, 075404 (2001). 

27 O. Giilseren, T. Yildirim, S. Ciraci and Q. Kihc, Phys. Rev. 
B 65, 1554XX (2002). 

28 D. Sanchez-Portal, E. Artacho, J.M. Soler, A. Rubio, and 
P. Ordejon, Phys. Rev. B 59, 12678 (1999). 

29 D.H. Robertson, DW. Brenner and J.W. Mintmire, Phys. 
Rev. B 45, 12592 (1992). 

30 G. G. Tibbets, J. Cryst. Growth 66, 632 (1984). 

31 K.N. Kudin, GE. Scuseria and B.I. Yakobson, Phys. Rev. 
B 64, 235406 (2001). 

32 H. Yorikawa and S. Muramatsu, Solid State Comm. 94, 
435 (1995). 

33 H. Yorikawa and S. Muramatsu, Phys. Rev. B 52, 2723 
(1995). 

34 C. J. Park, Y. H. Kim, K. J. Chang, Phys. Rev. B 60, 
10656 (1999). 



